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ABSTRACT 

It is generally considered that the neutron star cooling scenarios involving fast neu- 
trino emission, from a kaon or pion condensate, quark matter, or the direct Urea pro- 
cess, require the presence of baryon pairing in the central core of the star to control 
the strong neutrino emission and produce surface temperatures compatible with obser- 
vations. I show here that within the kaon condensate scenario pairing is not necessary 
if: 1) the equation of state is stiff enough for the star to have a thick crust in which 
sufficient friction can occur to heat the star and 2) a thin layer, of mass AM larger 
than ~ 10 -12 M Q , of light elements (H and He) is present at the stellar surface. The 
effect of the light elements is to increase the heat flow and thus produce a higher surface 
temperature. Both the occurrence of heating and the presence of H and/or He at the 
surface (deposited during the late post-supernova accretion) can possibly be confirmed 
or infirmed by future observations. 

Subject headings: stars: neutron dense matter stars: X-rays 

pulsars: general 
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1. INTRODUCTION 



2. OBSERVATIONAL DATA 



There are indications that some young pulsars may 
have surface temperatures lower than what is pre- 
dicted by the 'standard' model of neutron star cooling 
(see, e.g., Nomoto & Tsuruta 1986). Alternative mod- 
els with fast neutrino emission (see, e.g., Page 1994), 
however, predict temperatures much lower than ob- 
served. It had been proposed (Page & Baron 1990; 
Page & Applegate 1992) that baryon pairing (i.e., 
neutron superfluidity and/or proton superconductiv- 
ity, which strongly suppress the neutrino emission) in 
the densest regions of the neutron star core may be a 
way to keep the surface temperature high enough to 
be compatible with observations. Dissipation of the 
pulsar rotational energy into heat is another way to 
control the fast cooling, but detailed modeling of this 
effect (Shibazaki & Lamb 1988; Umeda et al 1993; 
Van Riper et al 1995) showed that it is by itself not 
sufficient. If the reported temperatures are close to 
the actual ones and if the actual temperatures are re- 
ally due to the cooling of the neutron star then pairing 
at very high density seemed to be a mandatory ingre- 
dient for the viability of the fast cooling scenarios. 

Recently, however, Potekhin, Chabrier & Yakovlev 
(1996) made an important point by noting that the 
presence of light elements (H and He) in the upper 
layers of the neutron star significantly affects the heat 
transport in these layers and results in a higher sur- 
face temperature when compared to the case of an 
iron envelope. They also showed that a layer of about 
10~ 12 Mq of accreted matter (possibly deposited dur- 
ing the post-supernova accretion phase) is sufficient to 
lead to quite different predictions for cooling models. 
I will show here that this hypothesis, when applied 
to a cooling model with fast neutrino emission from 
a kaon condensate (i.e., not as fast as the direct Urea 
process, Thorsson et al 1995), and including a rea- 
sonable amount of internal heating, leads to tempera- 
tures much higher than previously calculated and that 
'agreement with the data' can be obtained without the 
necessity of baryon pairing in the inner core where the 
fast neutrino emission is taking place, a result in total 
opposition to the previous models mentioned above. 

I summarize in § || the observational data on cool- 
ing neutron stars. My model is described in § || and 
the results are presented in § |]. The last section, ||, 
contains the discussion and conclusions. 



I summarize in this section the information presently 
available about estimates of and upper limits on the 
effective temperature of cooling neutron stars and I 
list in table [j] all relevant data. There are several more 
objects which I omit because the deduced upper limits 
on T e are still too high to be interesting. 

The output of cooling models is a luminosity L°°, 
i.e., an effective temperature T e °° related to L°° by 
L°° = AiraR 00 2 T e °° 4 (a being the Stefan-Boltzmann 
constant and R°° the neutron star's radius 'at infin- 
ity') and the observational problem is thus to deduce 
an accurate estimate of T e °° from the data. The quan- 
tities 'at infinity' are related to the local ones through 
the redshift factor (= 0.78 for my 1.4 M model) 
by: L°° = e^L, T°° = e*T, and = e^R. 
In the case of pulsars of goodness 1 to 3 (see ta- 
ble g) the rough upper limits are based on a com- 
parison with blackbody (BB) spectra. For the other 
five objects where there is some evidence in favor of 
the thermal origin of the detected photons (see, e.g., 
Ogelman 1995) the question of the spectrum used to 
obtain the temperature estimate is crucial. Model at- 
mospheres (e.g., Romani 1987; Shibanov et al 1992) 
show that the BB spectrum may be a very poor ap- 
proximation; in particular, magnetized hydrogen at- 
mosphere spectra at high temperatures, when the at- 
mosphere is wholly ionized, are much harder than 
BB spectra. At low temperature (< 10 6 K), how- 
ever, bound-bound and bound-free absorption may 
soften substantially the magnetized hydrogen spec- 
trum and push it toward or even below the BB spec- 
trum (Pavlov & Potekhin 1995), but no accurate cal- 
culations of magnetized spectra have yet been pub- 
lished in this regime. In the case of the Vela pulsar, 
the BB fits give a high temperature, T^ B ~ 1.6 • 10 6 
K, but require a neutron star radius between 3-4 km 
(Ogelman et al 1993) while fits with magnetized hy- 
drogen atmosphere spectra, which are expected to be 
quite accurate at such high temperature, give a tem- 
perature ~ 0.8 ■ 10 6 K, corresponding to a surface 
temperature Tjj ~ 10 6 K, for a neutron star radius 
of 10 km, R°° = 13.6 km (Page et al 1996). How- 
ever, if the luminosity deduced from the BB fits (with 
R°° ~ 3 - 4 km) is converted into an effective tem- 
perature T e °° for a neutron star with R°° — 13.6 km, 
then T e °° ~ TgP : this is the value I report in table [TJ. 
The neutron star 0002+6246 (Hailey & Craig 1995) 
is very similar to Vela and I calculate, in the same 
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way, T e °° from T^ B . In the other extreme, the cold 
Gcminga, presently available magnetized atmosphere 
spectra are too hard, and give very low temperatures 
which require an enormously large star or, alterna- 
tively, a very small distance (Meyer et al 1993), while 
BB fits (Halpern & Ruderman 1993) give higher tem- 
peratures which are compatible with emission from 
the whole surface of the neutron star at the known 
distance of 157 (+59, -34) pc (Caraveo et al 1996): 
I thus choose to use T^ B in this case, i.e., I hope 
that Tgg ~ r e °°. The two intermediate cases of PSR 
0656+14 and PSR 1055-52 are more delicate: it is dif- 
ficult to decide which of the BB or atmosphere spec- 
tra give temperature closer to the actual one. For 
PSR 0656+14, I report both values, TJ^ (Finley et 
al 1992) and TJf (Anderson et al 1993), the actual 
value of T e °° being probably somewhere in-between 
these, but for PSR 1055-52 only BB fits have been 
published (Ogelman & Finley 1993) so I use Tgg for 

± e ■ 

3. PHYSICS OF THE MODEL 

I will consider the fast neutrino cooling scenario 
induced by a kaon condensate (Brown et al 1988; 
Page & Baron 1990). The dense matter equation of 
state (EOS) is taken from Friedman & Pandharipande 
(1981), FP hereafter, and I assume that at densities 
above 

= 1.1 x 10 15 gm cm" 3 ~ 4 • p (1) 

the kaon condensate develops and induces a fast neu- 
trino emission 

^ = 10 24 (p/p ) 2/3 T$ ergs' 1 cm- 3 (2) 

(Thorsson et al 1995), where p is the matter density, 
po = 2.8 • 10 gm cm -3 the nuclear matter density, 
and Tg the temperature in units of 10 9 K. The critical 
neutron star mass for the appearance of a kaon con- 
densate is then 1.24 Mq . The modified Urea process 
and the similar bremstrahlung processes are active at 
all densities within the core but, at a temperature of 
10 s K, the kaon induced processes are more efficient 
by five orders of magnitude. Of upmost importance 
for the cooling is the presence and extend of baryon 
pairing in the core: the corresponding critical tem- 
peratures T c that I will use are plotted in figure [|. 
Notice that the assumed T c 's vanish within the con- 
densate region. In regions where pairing is present, 
the specific heat of the paired component is strongly 



reduced as well as all neutrino emission processes to 
which this component participates. 

Surface thermal emission is determined by the ef- 
fective temperature T e which is related to the interior 
temperature Tb through envelope models. T& is tradi- 
tionally taken at a density pj, = 10 10 gm cm -3 . The 
chemical composition of the envelope can strongly af- 
fect the Tfc — T e relationship: compared to models with 
pure iron composition, models with accreted matter 
give a higher T e for a given T&. I will use the recent 
calculations of Potekhin et al (1996) for both pure 
iron envelopes and envelopes with a layer of accreted 
matter, mostly H and He, of mass AM. For a 1.4 
Mq star with the FP EOS, the total mass of the en- 
velope, i.e., at densities below pb, is about 3 • lO -7 
Mq . I neglect here the effects of the magnetic field 
which are small (Page & Sarmiento 1996). 

Besides the neutrino and photon emissions which 
are energy sinks, rotational energy can be converted 
into heat by friction due to differential rotation of 
the inner crust neutron superfluid (Alpar et al 1984). 
I handle this effect in a simple way by writing the 
heating rate as 

ff « = J «' lo40 -(^)" S/2 < 3 > 

where t is the pulsar age, t = 300 yrs a typical 
spin-down time scale, and J44 the differential angu- 
lar momentum of the ff ictionally coupled crustal neu- 
tron superfluid in units of 10 44 g cm 2 rad s -1 . This 
expression assumes a standard spin-down rate from 
magnetic dipole radiation and is similar to the one 
used by Umeda et al (1993) and Shibazaki & Lamb 
(1988). This heating is distributed within the su- 
perfluid layers of the inner crust. I adopt the value 
J44 = 3.1 which corresponds to moderately strong 
heating, compatible with the size of the crust for the 
FP EOS. 

My treatment of other physical ingredients is de- 
scribed in previous works (Page & Baron 1990; Page 
& Applegate 1992; Page 1994) and all calculations 
were performed with a wholly general relativistic 
Henyey-type evolutionary code. 

4. RESULTS 

Figure || shows results for the cooling of a 1.4 
Mq neutron star without and with internal heating 
and various amounts of accreted mass AM, as well 
as for a 1.2 Mq star with heating and AM = 3 • 10 -8 
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Mq . As seen in figure [j], the 1.4 Mq star contains 
a kaon condensate, of mass 0.3 Mq , while the 1.2 
Mq star has a central density below p cr and thus fol- 
lows the 'standard' cooling scenario. During the neu- 
trino cooling era, i.e., at age up to about 10 4 - 10 5 
yrs when the neutrino luminosity is much higher than 
the surface photon luminosity, the nature of the en- 
velope material has no effect on the cooling rate but 
does affects the effective temperature: models with 
accreted envelope are warmer than models with iron 
envelope even though they all have the same interior 
temperature (considering separately models without 
heating and models with heating, the latter having of 
course a higher interior temperature than the form- 
ers). Even such a small amount of accreted mass as 
10~ 12 Mq has a very noticeable effect (Potekhin et al 
1996). When photon cooling takes over, the cooling is 
accelerated in the cases of accreted envelopes because 
of the increased surface emission. Finally, notice that 
at the latest ages the cooling is entirely controlled by 
the heating, independently of the previous history or 
envelope composition: the photon luminosity is sim- 
ply equal to the heating rate. 

Comparing these cooling curves, when both heat- 
ing and an accreted envelope are taken into account, 
with the plotted surface temperatures, one sees that 
the 'standard' cooling scenario, i.e., M < 1.24 Mq in 
my model, seems to be appropriate for PSR 1055- 
52 and maybe also PSR 0656+14 if the higher tem- 
perature estimate is considered for the latter. The 
temperature 'measurements' of Vela, NS 0002+6246, 
Geminga, and PSR 0656+14 if the lower T e value of 
the latter is considered, are below the standard cool- 
ing curve but can be fitted very well with the kaon 
condensate scenario in presence of heating and with 
an accreted envelope of mass larger than ~ 3 • 10~ 12 

Mq. 

5. DISCUSSION AND CONCLUSIONS 

The results of the previous section show that agree- 
ment between cooling models with fast neutrino emis- 
sion and the best presently available data does not 
necessarily require the presence of baryon pairing in 
the inner core of the neutron star: the superfluid- 
ity T c 's assumed in the calculations (figure |l|) vanish 
within the kaon condensate region. The three basic 
ingredients needed to obtain temperatures 'in agree- 
ment with the data' are: 1) fast neutrino emission as 
supplied by a kaon condensate, 2) moderately strong 



internal heating and 3) the presence of light elements 
in the envelope. 

There is no consensus on which, if any at all, fast 
neutrino emission mechanism is possible in neutron 
star cores but the kaon condensate is the least ef- 
ficient of all: a pion condensate or the direct Urea 
process (Thorsson et al 1995) would shift downward 
the cooling curves during the neutrino cooling era and 
require much higher heating rates. 

The internal heating by friction of the inner crust 
differentially rotating neutron supcrfluid depends on 
poorly known microscopic parameters (see, e.g., Van 
Riper et al 1995). A much stronger heating rate 
than assumed here is theoretically possible but would 
conflict with the temperature upper limits of PSR 
1929+10 and 0950+08 (see figure | and Van Riper 
et al 1995). Efficient heating, moreover, requires a 
thick enough crust, i.e., a rather stiff EOS, and there 
are arguments for such an EOS precisely in the case 
of the Vela pulsar from the analysis of its glitches 
(Link et al 1992) and from the effect of gravitational 
lensing on the amplitude of the pulses seen in its sur- 
face thermal emission (Page & Sarmiento 1996). An 
extreme softening of the EOS in presence of kaon con- 
densation was the result of the first models (Thors- 
son et al 1994) but the latest results of Fujii et al 
(1996) show that relativistic effects strongly reduce 
the softening of the EOS by the condensate: a stiff 
EOS even in the presence of a kaon condensate is thus 
not unlikely. Finally, in presence of heating the inde- 
pendence of the late cooling on the previous history 
and envelope structure could allow a determination of 
the microscopic parameters involved in the problem 
when detection of thermal emission from old pulsars 
will be achieved (Van Riper et al 1995). Actually, 
recent Hubble observations (Pavlov et al 1996a) may 
have detected the optical - UV Rayleigh-Jeans tail 
of the surface thermal emission of PSR 1929+10 and 
0950+08 and indicate surface temperatures of the or- 
der of 1 - 3 TO 5 K for the former and 6-8 TO 4 K for 
the latter. This estimated low temperature of PSR 
0950+08 is much below our simple model prediction 
but well above any cooling curve which does not in- 
clude any heating mechanism: if confirmed, it would 
show that heating is required but that the simplistic 
formula [j] is obviously much too naive. 

Due to the enormous uncertainty about the value 
of the pairing critical temperature T c at high density 
(see, e.g., Page 1994 for a discussion) it is a relief that 
the fast cooling scenarios, preferentially the kaon con- 
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densate one, do not necessarily require the presence 
of superfluidity. Takatsuka & Tamagaki (1995) re- 
cently emphasized that in presence of a kaon conden- 
sate both neutron and protons would pair in a 3 P2 
state, because of the high proton fraction, and nu- 
cleon pairing is then unlikely unless the neutron, or 
proton, effective mass is unrealistically high. The rel- 
ativistic effects described by Fujii et al (1996) imply 
precisely a strong reduction of the nucleon effective 
mass within the condensate region and disfavor pair- 
ing. However these effects also reduce significantly 
the softening of the EOS, implying thicker crust and 
hence higher heating rate, and the low nucleon effec- 
tive masses also reduce (Brown et al 1988). 

Finally, the main question arisen by the results 
presented above is the possible presence of light ele- 
ments at the surface of a neutron star. The late post- 
supernova accretion could deposit the needed 10~ 12 
Mq of light elements. The only indication I know 
of about the presence of hydrogen is in the case of 
the Vela pulsar when comparing spectral fits with BB 
spectra vs. magnetized hydrogen atmosphere spectra 
as discussed in § || which clearly show the inadequacy 
of blackbody-like spectra. Magnetized iron atmo- 
sphere spectra (Rajagopal et al 1996) are blackbody- 
like when used to fit ROSAT's PSPC low energy res- 
olution spectra: they would imply, for Vela, tempera- 
tures and stellar radii similar to the BB fits, i.e., radii 
at infinity R°° of the order of 3 - 4 km compared to 
14 km for magnetized hydrogen fits. Multiwcvclcngth 
studies (Pavlov et al 1996b) and the future X-ray mis- 
sions like SRG, AXAF, and XMM, will help resolve 
this issue. 

I have thus replaced the 'magic word' superfluid- 
ity by heating and accreted envelope: apparently one 
of them is needed to make the fast cooling scenarios 
viable. The presence of heating and/or of accreted 
matter will be confronted by observations in the near 
future while superfluidity has probably little chance 
of getting out of the theoretical realm. 

As a last comment I must emphasize that the 'clas- 
sification' of the observed neutron stars as follow- 
ing the slow ('standard') cooling scenario or some 
fast cooling scenario is model dependent. For ex- 
ample, the present results imply that Geminga and 
PSR 0656+14 (if the lower T e is used) belong to the 
fast cooling scenario but different assumptions about 
the extend of baryon pairing in the core can accom- 
modate them within the 'standard' cooling scenario 
(Page 1994). Moreover, Vela and PSR 0002+6246 



could also be accommodated within the 'standard' 
cooling scenario with a radically different assumption 
about the efficiency of the modified Urea rate (Schaab 
et al 1996). 

I thank A. Y. Potekhin for discussions as well as A. 
Y. P., G. Chabrier and D. G. Yakovlev for sending me 
their envelope models before publication. The idea of this 
work originated during the ECT* workshop 'Physics of 
Supernovae and Neutron Stars' (Trento, Italy, 3-14 June 
1996). This work was supported by a UNAM-DGAPA 
grant No. IN105495 and a Conacyt grant No. 2127P- 
E9507. 



Fig. 1. — Critical temperatures T c for neutron and 
proton pairing in the and 3 P2 channels used in 
this work (Page 1994). Indicated are also the cen- 
tral densities of the 1.2 Mq and 1.4 Mq neutron star 
models used and the crust-core and kaon condensate 
boundaries. 

Fig. 2. — Thin dotted line: Cooling curves for a 
1.2 Mq neutron star ('standard' cooling) with heat- 
ing and accreted envelope (AM = 3 ■ 10~ 8 ). Thick 
lines: cooling curves for a 1.4 Mq neutron star with 
a kaon condensate, various accreted masses and with 
or without heating. The data are from table 1. 
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Table 1 

Observational data on neutron star surface temperature 



NS 


Log Age 


Distance 


Log T e °° 


Goodness b 


Reference 




(yrs) 


(kpc) 


(K) 






0531+21 


3.10 


2. 


< 6.25 


3 


Becker & Ashenbach 1995 


0833-45 


4.05 


0.5 


5.87-5.91 


5 


Page et al 1996 


11 


11 




5.81 - 5.94 


5 


Ogelman et al 1993 c 


0002+6246 


~ 4 


- 3 


5.72 - 5.85 


4 


Hailey & Craig 1995 c 


1706-44 


4.24 


1.8 


< 6.1 


2 - 3 


Becker et al 1995 


2334+61 


4.61 


2.5 


< 6.1 


2 


Slane & Lloyd 1995 


1916+14 


4.95 


1.5 


< 6.0 


1 


Slane & Lloyd 1995 


0656+14 


5.04 


0.5 


5.93 - 5.97 


5 


Finley et al 1992 


11 


11 


0.3 


5.70 - 5.76 


J! 


Anderson et al 1993 


0740-28 


5.20 


1.9 


< 6.0 


1 


Slane & Lloyd 1995 


1822-09 


5.37 


1.0 


< 5.85 


1 


Slane & Lloyd 1995 


0630+178 


5.48 


0.16 


5.60 - 5.70 


5 


Halpern & Ruderman 1993 


1055-52 


5.73 


0.5 


5.74 - 5.82 


5 


Ogelman & Finley 1993 


1929+10 


6.49 


0.17 


< 5.5 


3 


Yancopoulos et al 1993 


0950+08 


7.23 


0.13 


< 5.1 


2 


Manning & Wilmore 1994 



a Listed are: neutron star name; log of age; distance; log of effective temperature; 'goodness' b of the data; 
reference. 

b Data 'goodness': 1: neutron star not detected, 2: neutron star detected at low count rate which precludes any 
serious analysis of the origin of the photons, 3: neutron star clearly detected but with evidence that the photons 
come mostly, or even exclusively, not from surface thermal emission, 4: neutron star clearly detected with some 
spectral evidence about the thermal origin of the photons and, finally, 5: neutron star clearly detected with good 
spectral evidence about the thermal origin of the photons. 

c The T e °° reported here is really an effective temperature obtained from the luminosity as discussed in section |^ 
and not the blackbody temperature T^ B reported by the authors. 
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